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We study discrete quantum dynamics where single evolution step consists of unitary system 
transformation followed by decoherence via coupling to an environment. Often non-Markovian 
memory effects are attributed to structured environments whereas here we take a more general 
approach within discrete setting. In addition of controlling the structure of the environment, we are 
interested in how local unitaries on the open system allow the appearance and control of memory 
effects. Our first simple qubit model, where local unitary is followed by dephasing, illustrates how 
memory effects arise despite of having no-structure in the environment the system is coupled with. 
We then elaborate this observation by constructing a model for open quantum walk where the unitary 
coin and transfer operation is augmented with dephasing of the coin. The results demonstrate 
that in the limit of strong dephasing within each evolution step, the combined coin-position open 
system always displays memory effects and their quantity is independent of the structure of the 
environment. Our construction makes possible an experimentally realizable open quantum walk 
with photons exhibiting non-Markovian features. 

PACS numbers: 03.65.Yz, 42.50.Lc 


I. INTRODUCTION 

In recent years, there has been a growing interest on defining and quantifying quantum memory effects in open 
system dynamics [iHini- These are based on a number of different approaches, ranging, e.g., from the concepts of 
information flow [2], to non-divisibility and mutual information [5]. Based on these theoretical developments, 

experimental realizations for detecting and controlling of non-Markovianity has become feasible mm and also 
a number of proposals to exploit memory effects, e.g., to quantum information tasks has been recently proposed 
[niiaiii]. The previous studies on non-Markovian quantum dynamics mostly focus on continuous coupling between 
the open system and its environment. However, there also exists other possibilities to study memory effects, for 
example discrete dynamics which we consider here. In this case, the open system couples in stepwise manner in time 
to its environment. One can also envisage a possibility of having unitary transformations changing the state of the 
open systems between non-unitary evolution steps caused by the environment. We are interested in what is the state 
of the open system after each step which consists consecutive unitary and non-unitary parts. Thereby, the theoretical 
challenge is to construct a discrete dynamical map for the stepwise evolution of the system of interest and this map 
should contain information about the local unitary within the system of interest and the properties of the environment 
it is interacting with. 

Quantum walks provide a promising base to combine the study of discrete dynamics and memory effects. Generally, 
in quantum walks the system can evolve discretely or continuously nMii] and a study on the relation between the two 
cases can be found in m- Quantum walks have proven to be important in such diverse fields as quantum information 
processing m, complex networks [20l[2T] and the physics of topological phases [22l[23j. Moreover, the dynamics of 
quantum walks show a very broad range of different dynamical behavior from ballistic spreading to localization m- 
132] . Introducing noise to the quantum walk, the effects of the transition from unitary to non-unitary dynamics and 
the classical limit can be studied [mnaiMiisMi- Recently, quantum walks have been implemented experimentally 
for example using trapped ions [47l |48j , atoms in optical lattices |49] , linear optics [23l [37] , optical fibers [28l |50| and 
waveguides pSlISTl 15^ . 

In this article we study non-Markovian discrete dynamics of discrete quantum walk in a line. For this purpose, we 
introduce a new type open quantum walk where the dynamics is characterized by a non-divisible discrete dynamical 
map. It is worth pointing out that, to the best of our knowledge, all previous approaches to decoherent quantum walks 
can be characterized by using discrete quantum dynamical semigroup. We go beyond dynamics which is describable 
by semigroup and are interested in how to induce and control memory effects for quantum walks. The article is 


kimmo.luoma@tu-dresden.de 




2 


structured in the following way. We first introduce the concept of discrete dynamical map and also formulate the 
suitable quantifier for memory effects. To understand better the origin of memory effects in the considered models, we 
then review the standard dephasing model for a qubit with non-Markovian dynamics. We then study discrete open 
qubit dynamics with unitary control operation and show how the addition of the control transforms the dynamics of 
the open system inherently non-Markovian. The elaborate on the insights obtained from the simple qubit model, we 
construct a discrete quantum walk where we identify the coin operator as the local control operation and proceed with 
memory effects in 1-d discrete open quantum walk. Throughout this paper the interaction between the open system 
and its environment is of pure dephasing type and in the quantum walk model the coin is coupled to an environment. 
The formulation follows the path and emphasis the experimental realization of the models with photons. 


II. DISCRETE DYNAMICAL MAP 

Quantum dynamics is generated by quantum dynamical maps, a family of completely positive and trace preserving 
(CPT) maps , n > 0, such that $0 = I • If the state of the quantum system, described by a density operator or 
matrix, is initially po then = d)„(/9o) defines the evolution of the state. Discrete dynamics for a quantum system 
emerge if the parameter n indexing the family of CPT maps takes discrete values only, eg. n G N. We assume 
that 0 G N. 

CPT maps suitable for studying discrete dynamics can be constructed by using the Stinespring’s dilation theorem 
[55] . It states that for every CPT map it is possible to assign an unitary evolution in enlarged Hilbert space. Total 
Hilbert space is thus H = 'Hs®'Hs and the map is generated as pn = tr£:{[/"poC)x(17l)"}. Here U is unitary operator 
acting on total Hilbert space and x fixed initial state on the auxiliary Hilbert space. This construction guarantees 
that is CPT for each n. Physical interpretation for this construction in the context of open quantum systems is 
that T-Ls is a Hilbert space for external environment to which the system is coupled unitarily. 

Quantum dynamical maps can be classified in the following way by their divisibility properties. If the dynamical 
map satisfies the following decomposition law 

® ( 1 ) 

for all n G N, then the dynamical map is Markovian and forms a discrete dynamical semigroup. The dynamical map 
is called CP-divisible if 


^n+m — IPn-t-m,m ® ( 2 ) 

for n, TO G N and where Wn+m,m is completely positive and trace preserving. If the dynamical map can be composed 
as 


^n+m —I^,m (3) 

where Vn,m is positive and trace preserving then the dynamical map is called P-divisible [54] . 

It is clear that the divisibility property of the dynamical map goes beyond the concept of semigroup. However, 
the various definitions of quantum non-Markovianity are still under active discussion m- In the next section we will 
introduce the measure used in this work to quantify non-Markovianity of the discrete dynamical map. 


III. MEASURE FOR NON-MARKOVIANITY 

Following |5], we use a measure based on distinguishability of quantum states, quantified by the trace distance 
d{pi,P 2 ) = Itrjpi — P 2 \- Trace distance is contractive under positive and trace preserving maps [55], hence also under 
completely positive and trace preserving maps. Thus, the evolution of trace distance for fixed initial pair under P- 
divisible map is contractive, which means that Dp^ p^{m) > Dp^ p^(n) for all to < n where Dp^ p^(n) = d{^nPi, ^nP 2 )- 
If we find that the trace distance increases, Dp^^p^(n) > Dp-^^p^(rn) for some n > m, then we know that the dynamical 
map is not P-divisible and we say that the dynamics is non-Markovian. 

To construct a measure for discrete dynamics, we define the increment of the trace distance evolution 


Au 2 (n) =Dp^^p^{n) - Dp^^p^{n -l),n> 1, 
Ai,2(0) =0, 




3 


and the measure for non-Markovianity is defined in terms of the increment as 

A/’($)(n) = max Ai, 2 (n) 

Pl',P2 

nGS={nGN| Ai, 2('JT')>0} 


( 4 ) 


It can be shown that it is sufficient to make the maximization over orthogonal pairs of states [56] and it also has a local 
representation HU. The measure has physical interpretation in terms of information flow between the system and 
the environment. When Ai^ 2 (?^) < 0 information flows away from the system to the environment and the ability to 
distinguish the two states decreases. When Ai^ 2 {n) > 0 there is a backflow of information from the environment to the 
system which improves the distinguishability. In general, a lower bound for non-Markovianity is quite straightforward 
to obtain using only a small number of initial states. It is also worth noting that there exists evolutions which are not 
P-divisible but the trace distance between evolving states might still decrease at all points of time. In this work we 
do not present the results for full optimization over all the initial states but instead focus on specific pairs of initial 
states which allow to witness the non-Markovianity of the dynamics. 


IV. DISCRETE QUBIT DYNAMICS 
A. Dephasing 


Genuine quantum effect on open quantum system dynamics is the loss of quantum coherences without energy 
exchange between the system and the environment. This effect is called pure dephasing. Our motivation to study this 
effect on qubit dynamics is due to possibility of experimental implementation using optical elements |12j . Coupling 
between the system and the environment is given by 

Ust = f doj ( 5 ) 

'' v=L,R 

where \i') labels the qubit degree of freedom and jw) corresponds to the environmental degree of freedom. With an 
optical implementation in mind, this type of unitary dynamics describes the interaction of a photon with a birefringent 
medium, eg. quartz m- Different basis states \v) correspond to the different polarization states, rii, is the polarization 
dependent index of refraction, ui is the frequency of the photon and 5t corresponds to the thickness of the quartz 
plate. 

For a fixed product initial state g = p{d)®x this coupling generates pure dephasing dynamics, given by a dynamical 
map <i>™p(0) = tr£ expressed in \v) basis as 

( \iy){v\ ho Vi/, 

J \L){R\ ^n{St)\L){R\, (6) 

[ |i?)(L| ^R*mR){L\ 

The map is determined by the the function k : K_|_ i—>■ C, K{St) = J du;e*^"“'^*|x(w)P, where An = — rin. 

Throughout this work we take the environment initial population distribution |x(a;)p (spectrum) to be constructed 
of two Gaussians with widths tr, amplitudes and where A G [0,1], central frequencies of the peaks pi,p 2 
and peak separation 5uj = p 2 — Pi, so that 

=1^ . (7) 

This choice is experimentally motivated |12j and with parameters A and a it is possible to control the structure of 
the environment. 

It can be shown that the dynamics over a period 6t can be Markovian or non-Markovian depending on the properties 
of the population distribution |x(a;)p of the initial state of the environment. Dynamics as a function oi St = t is 
plotted in Fig.[^ Trace distance dynamics for the optimal pair of initial states is given by Dp^ p^{5t) = |K((5t)|. This 
type of dynamics is analyzed in HH. 
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FIG. 1. Standard dephasing model for a qubit. On the left is plotted the spectrum of the environment, it consists of Gaussian 
peaks. On the right is plotted the trace distance dynamics for maximizing initial state pair. As one can see, the single Gaussian 
peak gives Markovian dynamics and more structured environment gives rise to quantum memory effects. Here = 5uj/{2n). 


B. Dephasing with local control 

We consider now discrete dynamics where at each step in the dilatation space we have an action of a local control 
unitary operation m followed by dephasing. The single step unitary in the total space is then 


V = Usf{Cr,®h). ( 8 ) 

In this work we consider only the following local control unitaries Crj 

C, =y^{\L){L\ - \R)m + ^^^{\L){R\ + |i?)(L|). (9) 

These correspond to biased beam splitter transformations and = 5 being the Hadamard transformation also known 
as the balanced beam splitter transformation. The local unitary operator transforms the polarization basis into a new 
basis which is not generally simultaneously diagonalizable with the decoherence basis. This proves to be crucial for 
the non-Markovianity of the quantum dynamics as we will show. The reduced dynamics generated by is given by 
which is defined as 


p(n) = ^^p=tvs {V'^{p®x){y'^T] ■ (10) 

The state of the qubit can be expressed as p = 1 (I 2 + r ■ a), where r = {ri,r 2 ,ry^ is the Bloch vector and 
(7 = (ui, (72, where ai are the usual Pauli matrices. For arbitrary values of 77 the dynamics is most conveniently 
given in terms of the Bloch vector. The dynamical map can be written as 


rim 




( 11 ) 

( 12 ) 
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The matrix M(uj) is given by 

( —f3 cos AnStoj — sm{AnStuj) a cos AnStuj \ 

PsinAnStu; —cosAnStuj —a sin Arn^to; , (13) 

a 0 ^ / 

where a = 2 a / (1 — 77 ) 77 , /3 = 277 — 1 and An = ul — n^. Matrix M(w) is periodic, eg. M{u} + 11 ) = M{uj), where 
tl = 5^77 • Numerical integration of M(a;)’” must be done carefully since integrands are highly oscillatory. However, 
reliable numerical results can be obtained. In the strong dephasing limit, a and for some particular values of rj 
it is possible to obtain analytical results. 

Figures and present the results for the values rj = 0,0.5,1.0 in the case of weak and intermediate dephas¬ 
ing strength. It is worth noting that for all of the cases we have A = 0, i.e., the environmental spectrum is flat 
corresponding to Markovian dynamics without local control. 

For 77 = 1, the local control is Ci = cjz and we have 

{n,u:\iUst ■ {az®l£)r\L,uj) 

{v,uj\{Ust ■ {oz®l£)r\R.uj) 

This leads to the following dynamical map 

\L){R\ K{m5t){-ir\L){Rl (16) 

[ \R){L\ ^ K*{m6t){-1)'^\R){L\. 

Compared to the uncontrolled case, the sign of the coherences is flipped. This does not affect the memory effects since 
the dephasing process moves the states towards the z-axis in the xy-plane of the Bloch sphere and the sign change of 
the coherences keeps the distance of the state from the z-axis constant. The dynamics is displayed in panels c) and 
d) of Fig. Hand Fig. § 

For 77 = 0 the local control is Gj. and we have 

{n,u:\{Ust ■ (a,®If))“/',a;) =(i _ ( 18 ) 

This gives the dynamical map 


(14) 

(15) 


(19) 

{\v){v\ ^\v'){u'iyv^v', 

<^'^^+^:{\L){R\ ^K{5tr\R){L\, ( 20 ) 

[ \R){L\ ^ K{5t)\L){R\. 

In this case, the information flow between the system and the environment is maximal in the sense that the local 
control is able to completely eliminate the effect of the environment after even number of steps, see panels a) and d) 
in Fig. and Fig. 


C. Strong dephasing limit 


Let us assume that we have a “flat” population distribution for the initial environmental state. What we mean by 
this is that, \x{oj)\^ stays almost constant for w e [w' — f2/2,a;' + f2/2], where w' G [0, 00 ). 

Now, this means that the effects of the environment in this limit are generic in the sense, that the global structure 
of |x(w)p does not play a role. This type of environment is usually called ’Markovian’. It is intuitively clear that 
in this type of situation it is enough to integrate only over a single period of length Q. in Eq. (Ill to a very good 
approximation. In terms of equations this means 


1 

r{m) =— / dujM{uj)'^f. (21) 

Jo 
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FIG. 2. Weak dephasing interaction. Panel a) r] = Q,h) rj = 0.5, c) 77 = 1.0 and we have plotted the time evolution of non-zero 
Bloch vector components. Panel d) plots the measure for non-Markovianity for fixed pair of initial states for the used three 
values 77 . Initial states in all figures are ri = ^(1,0,1)^ and r 2 = —f\ expressed in terms of the Bloch vector. Parameters 
of the environment are ^ = 0 and Su) = 9a. Parameters for the interaction are An = 0.009 and St = 0.014 ^^^^ . This gives 
§ « 643 » 1. 


Physical idea behind this approximation is that there is no distinguished frequency of the environment because of the 
flatness of the spectrum. Then all the frequencies in the interval of length 17 “mix” the state of the reduced system 
with equal weights, hence it is sufficient to consider only one of these intervals. 

In general we have to validate this approximation numerically, but for the special case where = 5 and ^ = 0 we 
can obtain analytical expressions for dynamical map in the strong dephasing limit, which we will do next. It turns 
out that we will need the condition a ^ fi for the analytical calculation but the numerical data shows that already 
for intermediate dephasing, ct ~ 17, strong dephasing approximation (21) works quite well, see Fig. Dynamics in 
the intermediate and strong dephasing regime for 77 = ^ is plotted in panel b) of Figs. and 




We 


We now give a more detailed proof of Eq. (21) for 77 = ^ and A = 0. Let 77 = ^ and |x(<^)P = 

assume that 17 <C cr, which means that the spectral distribution varies in much larger scale than the integral kernel 
M{uj), which takes the following form 


( 0 — sin(A 7 r< 5 ta;) cos{AnStoj) \ 

0 — coslAnStuj) — sm{An6tuj) j . (22) 

10 0 ) 

We decompose the interval K+ = [ 0 ,oo) as UfeGN\{o} where 17*, = [17 • (fc — 1 ),D • k). Since du;|x(a;)p is 

continuous on closed interval Dj, and differentiable on 17/^, mean value theorem states that we can always find oj^ € 
such that 


|x(a;fc)pf7= [ dw|x(w)p. 
J 


( 23 ) 
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FIG. 3. Intermediate dephasing interaction. Panel a) 77 = 0, b) 77 = 0.5, c) 77 = 1.0 and we have plotted the time evolution of 
non-zero Bloch vector components. Panel d) plots the measure for non-Markovianity for fixed pair of initial states for the used 
three values 77 . Initial states in all figures are ri = r 2 = —ri expressed in terms of Bloch vector. Parameters 

of the environment are ^ = 0 and Suj = 9a. Parameters for the interaction are An = 0.009 and St = 2 ■ This gives 

B = 4.5 ^ 1. 

(T 


As an approximation we choose the midpoint of each interval, cok = ^ ■ k + then 


(24) 

/cGN\{0} 


where 793 ( 11 ,( 7 ) is Jacobi theta function. 

Condition <C cr allows to do the following approximation 


du;|x(^)|"M(a;)™«^|x(wfc)|M dccM(u;)'" = ^ |x(a;fe)|2M(m) 

k k k 

= i^3(7r(i - /r/fi), e-2-^'^'/*^^)M(m), (25) 


where M{m) = dojMlujy 


and we used Eq. (241 in the last step. 


When we take the strong dephasing limit, 

. Next we will 


a —>■ 00 and use the property limq_,,o 793 ( 77 , g) = 1 of the Jacobi theta function we obtain Eq. (21) 
construct the dynamical map in this special case. 















D. Dynamical map in strong dephasing limit for A — Q and V ~ ^ 

In the strong dephasing limit for single Gaussian spectral distribution and = 5 we obtain the following analytical 


form for the dynamical map ^ /o 


Ao —I, 


^0 0 r 
A2 I 0 1 0 
,0 0 0, 


'0 0 o> 

Ai = I 0 0 0 I , 
.1 0 Oy 


1 


'I 0 r 


A3 =- 10 1 0 
,0 0 1, 



1 (Xfji—2 0 



1 (^m — 2 0 (^m—2\ 

A2m — j 

0 hm-i 0 

1 , TO > 2 , 

A2m-1 — j 

0 bm-2 0 j 


\^m—2 0 O.m—2/ 



3 0 ^m—2/ 


, TO > 3. 


(26) 


where 


C(i) 

Qk + 

cik =h = 0 , A: < 0 , 


Cik)= ' 


k + l\k 

C{k) is called a Catalan number, and bk have the following limiting behavior 


k€N. 


(27) 

(28) 


lim bk ={V2 — 1 ), 

fc—>-oo 


lim Ufe =1 - j=. 

fc —>00 y 2 


Thus the dynamical map takes the following form in the limit of inhnite number of steps 


(29) 


b 

1 

V2 

0 

^ V2 

■^m—>-oo — 1 

0 

v^- 1 

0 

1’ 

1 

vZ 

0 



(30) 


E. Discussion 

The results above show that in a discrete dephasing model for a qubit, the addition of local unitary can induce 
non-Markovian dynamics even for flat ’’Markovian” spectral structure. Local unitary operation can be seen as a 
periodic control, that dynamically decouples the open system from the environment giving rise to a partial revivals 
of the populations and coherences in the open system dynamics. Periodicity of this control operation can be seen 
from the discreteness of the dynamics, i.e., we ’’watch” the system only at the integer multiples of the control period. 
It is also worth noting that the local unitary changes the open system state while the earlier created correlations 
between the system and environment still persist. The local change in the system state allows then the existing 
system-environment correlations to be converted back to the increased distinguishability of the system states and 
backflow of information. 

Another effect of the local control, in the case that it does not commute with the dephasing basis, e.g. when 
Co = o'x, is that it transforms the open system dynamics from pure dephasing to dissipative. For the case when local 
control is Ci = < 7 z, which commutes with the dephasing operator Ugt, the dynamics is pure dephasing type and the 
action of the local control does not have effect on the non-Markovianity of the dynamics, i.e., non-Markovian dynamics 
can emerge from the spectral structure only. We also show that the non-Markovianity induced by the local control 
is generic in the sense that the structure of the environment spectrum does not play a role in the strong dephasing 
regime. In the special case of Hadamard control ry = 0.5 and A = 0 we were able to derive analytical expression for 
the dynamical map. 

In the following section we study more complicated situation with one dimensional discrete quantum walk. Our 
findings will be better understood with the help of the physical intuition gained from the present section. 
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FIG. 4. Error of the strong dephasing approximation for few values rj. Error is measured by — A„), eg. the trace 

distance between the exact and approximate (A„) dynamical map. For this calculation we have chosen parameters of 

the environment as A = 0 and Suj = 9cr. For interaction we have An = 0.009 and the solid lines correspond St = 0.02 ^^^^ 
(weak dephasing) and dashed lines correspond to St = 1.03 (intermediate dephsing). These corresponds to the following 
approximate values Q./a k, {450, 9} of the quantity measuring the strength of the dephasing. 


V. OPEN QUANTUM WALK 
A. Quantum walk 


Quantum walks are either continuous or discrete time unitary protocols that evolve a quantum state on a Hilbert 
space that is constructed from the underlying graph where the walk takes place. In this section we will mostly adapt 
to the notation of Ref. . 

In this work we limit the discussion to a discrete quantum walks on a line. Hilbert space for the walk is T-Lw = 
"He 0 "Hp = 0 €^(Z). It consists of the coin- and the position space. Unitary operator W evolving the state of the 

walker over a single step is the following 

W =i\L){L\ 0 5+ \R){R\ 0 5t)(Cp 0 Ip), (31) 

where S = J2x 1^ ~ l)(^l = Ci , see Eq. We choose to focus only on Hadamard walks, meaning that the 

coin unitary is the Hadamard matrix. The unitary operator W can be diagonalized if we move into a quasi-momentum 
picture by Fourier transform \k) — k G [—7r,7r). k is called quasi-momentum since it is periodic. Note 

that the Fourier- or quasi-momentum basis is not normalizable, but nevertheless very useful when used carefully. The 
inverse transformation is defined as \y) = ^ J^^dke~^^^\k). 

In this work we always initialize the position of the walker to the origin. In the quasi-momentum picture the unitary 
operator W acts on an arbitrary state initialized from origin, |(/)o) = |0) 0 |0), as 

dA- 

— (Mfe)'"|(/.)0|fc), 


( 32 ) 
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where Mk is the following 2x2 matrix 


Mk 


1 

71 



Eigenvalues of Mk are {e *'^'“,—6®'^'=}, where i^k is defined by 


sinfc =-\/2sinj^fe. 


(33) 


(34) 


Using the quasi-momentum representation for solving the dynamics and then transforming back to the position 
representation, we obtain the following expression for a general initial state starting from the origin, I'lpo) = (cl|E) -|- 
Cij|i?)) 0 |0), and evolving m steps 


\^m) = = E 


x——m 


+ cnA^ix))\L) + {clB^{x) + chB)?(x)) \R) 


0 |x). 


(35) 


Analytical expressions for the coefficient functions A™(x), A^(x), i3£®(x), B^{x), are obtained most easily from the 
quasi-momentum picture. They are 


where 


Al^{x) =—^ - [a^{x) + r{x% 

B^{x) = -7L^Z-.^ [/3"®(c,)+y-(x)], 


1 (_l)m+x 

A^{x) =—^ - [r{x) - 7 “(x)], 

1 (_l)m+x 

B^{x) =-' ^- A c^^ix) + /3”®(x)], 


a”®(x) 

rix) 

7™(x) 



dk 

2 tt 


i(kx—muk) 
^ •) 


dk cos k 

27r 1 -I- cos k 

dk sin k 

27r -^/l -I- cos k 


^i{kx—mvk 

^i(kx-mvk 


(36) 

(37) 

(38) 

(39) 

(40) 


Using the method of stationary phase, asymptotic expressions for Eqs. (36) and (37) may be obtained. 

After m number of steps the walker has non-zero probability to be found in positions —m, —m + 2, • • • , m — 2, to. 
From this follows that after even number of steps the walker has non-zero probability only at even vertices, similarly 
for odd number of steps. This property is sometimes called modularity. 


B. Open quantum walk 

We extend the the Hilbert space to take into account the effect of an external environment. The extended Hilbert 
space is 77 = He ®Hp ® He = Hs ® He and now our system consists of the coin and position degrees of freedom. 

We couple the coin degrees of freedom to the external environment and the coupling unitary is given by the trivial 
extension of Eq. ([^ 

Ust=[duj ^ (8)Ip (g) |a;)(a;|. (41) 

For simplicity we consider only homogeneous case, meaning that the coupling operator acts identically on each vertex 
of the graph. The walk operator is extended to act on the enlarged Hilbert space trivially, W = W®le. 

The dynamical map for the open quantum walk is constructed as usual 

4>)^(po) =tr£ {{Ust ■ Wrp^®x{{Ust ■ W)^r} , (42) 

where x is again initial state for the environment. Remarkably, an analytical expression for the dynamical map can 
be found. 
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1. Quantum dynamical map 


Matrix elements of the dynamical map are [^^]a'a,x'x-,TT',yy' = {<^\x'\^^{\a,x){T,y\)\T',y'). We use the following 
shorthand notation [<i)^]o-/cr,x';rT',y' = [^m]cr'iy,x'o-,TT', 0 y' ■ Matrix elements can be written explicitly by using the 
definition of the quantum walk 


= /dwIxHp ^ 

,fci={i,fi} 

In,---,li={L,R} 

^ (^Cln,ln — 1 * * ’ |^n)(^n|^ ) 

{x'\K,^---Kk,\0){0\Kl---Kl\y'), (43) 


where = S, Kji = and Ca,a' = {<^\Ch\<^')- 

Each non-zero inner product {x'\Kk^ ■ ■ •it'fcjO) in Eq. (43) corresponds to a path 0 
Each path 0 


of Eq. (43) as 


x' in a n-level binary tree. 
X has N'l{x) = and N'^{x) = left and right turns. This allows to write the exponential part 


_ iujStAn{y'-x') 


(44) 


where An = ul — ur, since left and right turns each contribute to the sum in total Nr{x) and NR{n) times the term 
riL and ur, respectively. Then by suppressing the explicit matrix products we can write the matrix element as 




a' (T,x' ;tt' ,y' 


= / da;|x(a;)|2e*‘"'^‘^”(^'-"')(a',x'|W”|a,0)(r,0|(W^)"|r',y')- 


(45) 


This shows, that the coupling of the coin degrees of freedom to the external environment has an effect only on the 
coherences between different sites. It does not have any effect on the position distribution or to the propagation speed. 


2. Strong dephasing limit 

In the strong dephasing limit the state of the quantum system is block diagonal since the coherences between any 
two different sites in the lattice are destroyed. This means that the state of the quantum walker can be written as 

Pm= 0 P^{x), (46) 

x^[—m,rn\ 


where p^{x) = Yha (j'={l R) Eigenvalues of a self-adjoint 2x2 matrix are A± = 

^ (^a + b± Y^(a(^^cp~T4j5|2^, 

Pure initial state IV'o) = cl|L, 0) -I- Ci^|i?,0) evolves in m-steps in the strong dephasing limit into 4>)^(|i/)o)(V'oD- 
Expression for block x is 


{x\^Zi\^o){^o\)\x) = (a;|IT’"|V'o)(V'o|(IEl')’”|a;) = p^oix). 


(47) 


Characterization of non-Markovianity becomes now more feasible, since we need to diagonalize at the mth step 2 to-|- I 
2x2 self-adjoint matrices, instead of valued matrix. Analytical expression by using Eqs. (35),(36),(37) could be 
obtained, but in this work it we do not analyze it further. 


3. Non-Markovianity of the open guantum walk 

We choose the initial state pair for the walk to be |T)|0) and |i?)|0) for simplicity. Parameters of the environment 
are Sco = 9a, An = 0.009, and we vary the value of parameter A which controls the structure of the environment. 
Non-Markovianity as a function of the dimensionless parameter StVlAn is plotted in Fig. 5 where SI = 

We can see that the measure for non-Markovianity is periodic when we increase 5t iraich gives the duration of 
the dephasing step. Remarkably, for a specific values of 5t the structure of the environment is irrelevant, i.e., for 
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FIG. 5. Non-Markovianity of the Quantum walk after ten time steps, A/io- Parameter A controls the structure of the 
environment spectrum. See text for other parameter values. Non-Markovianity is plotted as a function of the dimensionless 
interaction time of each dephasing unit. As StflAn increases the dephasing effect of each Ust becomes stronger and the solid 
black line corresponds to the strong dephasing limit. 


specific duration of dephasing, the quantity of memory effects remain the same irrespective of the structure of the 
environment. The periodicity of the measure is related to the separation of the peaks of the environment spectrum 
in the sense that A/’|a=o = J^lA^a when St = m ^ . It is also very interesting to note that for strong dephasing 
limit, the value of non-Markovianity approaches a constant value which is again independent of the structure of the 
environment which resembles the simple discrete qubit dynamics case. However, before the strong dephasing limit 
when St has small or intermediate value, non-Markovianity has two sources: the local unitary and the environmental 
structure. For A = 0, the non-Markovianity originates from the coin flipping only (local unitary) while for other 
presented values of A, the structure of the environment gives additional contribution to the quantity of the memory 
effects. 


VI. SUMMARY 


We have analyzed in detail two different models for discrete open quantum system dynamics: a simple qubit 
dynamics with dephasing and more elaborate full open quantum walk model. The results show that the addition 
of local unitary operation can dramatically change the nature of the open system dynamics and in particular the 
appearance of memory effects. We have discussed how memory effects generated within our models can become 
generic in the sense that the structure of the environment spectrum does not play a role - the amount of memory 
effects is independent of the form of the spectra - and non-Markovianity is solely induced by local unitary operation 
in these cases. The considered models and observed phenomena give rise to approximation schemes which simplify 
the description of the dynamics considerably that may prove useful for other contexts too. We have also discussed in 
detail for the first time a model for open quantum walk where memory effects can be introduced and controlled in an 
experimentally relevant way. The presented open quantum walk is experimentally realizable using linear optics. 
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